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Abstract 

The massive Schwinger model is studied, using a density matrix renormahza- 
tion group approach to the staggered lattice Hamiltonian version of the model. 
Lattice sizes up to 256 sites are calculated, and the estimates in the contin- 
uum limit are almost two orders of magnitude more accurate than previous 
calculations. Coleman's picture of 'half-asymptotic' particles at background 
field 9 = TT is confirmed. The predicted phase transition at finite fermion mass 
(m/g) is accurately located, and demonstrated to belong in the 2D Ising uni- 
versality class. 

PACS Indices: 12.20.-m, 11.15.Ha 
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The Schwinger model or quantum electrodynamics in one space and one time dimen- 
sion, exhibits many analogies with QCD, including confinement, chiral symmetry breaking, 
charge shielding, and a topological ^-vacuum IQ-Q. It is a common test bed for the trial 
of new techniques for the study of QCD: for instance, several authors have discussed new 
methods of treating lattice fermions using the Schwinger model as an example ||^. 

Our purpose in this paper is twofold. First, we aim to explore the physics of this model 
when an external 'background' electric field is applied, as discussed long ago in a beautiful 
paper by Coleman [Q. Secondly, we wish to demonstrate the application of density matrix 
renormalization group (DMRG) methods [^,[^ to a model of this sort, with long-range, non- 
local Coulomb interactions. The DMRG approach has been used with great success for 
lattice spin models and lattice electron models such as the Hubbard model; here we apply 
it to a lattice gauge model. 

Choosing a time-like axial gauge Aq = 0, the Schwinger model Hamiltonian becomes 

H = j dx(^-ii}-i^{di+igAi)^p + mi)i) + ]^E^^ . (1) 

where if) is a. 2-component spinor field, since there is no spin in one space dimension. The 
coupling g has dimensions of mass, so the theory is super-renormalizable. Using g as the 
scale of energy, the physical properties of the model are then functions of the dimensionless 
ratio m/g. Gauss' law becomes 

d^E = -d^A' = gf = (7V^7V. (2) 
where E is the 1-component electric field. This can be integrated to give 

E = gjdxf{x) + E, (3) 

showing that E is not an independent field, but can be determined in terms of the charge 
density j°(a;), up to the constant of integration F, which corresponds to a "background 
field", as discussed by Coleman [||. We can think of the background field as created by 
condenser plates at either end of our one-dimensional universe. 

If |F| > g/2, charged pairs will be produced, and separate to infinity, until the field is 
reduced within the range \F\ < g/2, thus lowering the electrostatic energy per unit length. 
Thus physics is periodic in F with period g, and it is convenient to define an angle 6 by 
6 = 2TtF/ g. Then we can always choose 9 to lie in the interval [0, 27r]. 

In the weak-coupling limit m/g — *• cxo, the vacuum contains no fermionic excitations, and 
the vacuum energy density eo is given purely by the electrostatic energy term (we ignore the 
energy of the Dirac sea) eo = F'^/2 {6 < tt) or {g - Ff/2 (tt <e <2'k). 

Thus there is a discontinuity in the slope of the energy density, corresponding to a first- 
order phase transition, at ^ = vr. In the strong-coupling limit m/g = 0, on the other hand, 
chiral invariance demands that the vacuum energy density remains constant as a function 
of 6' [0. Thus we expect a first-order transition at ^ = vr for large m/g, which terminates at 
a second order critical point at some finite {m/g)c- The two-fold vacuum degeneracy would 
naturally lead one to expect that the transition will lie in the Ising universality class. 

Normally, charge is confined in the model: there is a 'string' of constant electric field (or 
fiux) connecting any pair of opposite charges @J^. But Coleman points out that in the 
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very special case 6 = ir, the peculiar phenomenon of "half-asymptotic" particles arises. In 
the weak-coupling limit, one can envisage the state shown in Fig. The electric field energy 
density is the same in between each pair of particles, and they can therefore move freely, as 
long as they maintain the same ordering, i.e. no pair of fermions interchanges positions. 

The one-dimensional fermionic theory can be mapped into an equivalent Bose form . 
At m/g = 0, the theory is exactly solvable, and reduces to a theory of free, massive bosons, 
with mass M = gl\pji., independent of the background field Q. At large mj g and ^ = vr, 
where the spontaneous symmetry breaking occurs, the half-asymptotic particles are found 
to correspond to topological solitons, or "kinks" and "antikinks" , in the Bose language |Q . 

There have been very few attempts to verify these predictions numerically, that we are 
aware of. Hamer, Kogut, Crewther and Mazzolini used finite-lattice techniques to address 
the problem, and located the phase transition at 6' = vr to lie at [mj g\ = 0.325(20), with 
a correlation length index u = 0.9(1). Schiller and Ranft |Tl| used Monte Carlo techniques 



to locate the phase transition at {m/g)c = 0.31(1). For zero background field, on the other 
hand, many different numerical methods have been applied to the Schwinger model; for a 
recent review, see Sriganesh et al. [p!0[|. 



We employ the Kogut-Susskind |12| Hamiltonian lattice formulation of the Schwinger 



model, following Banks et al. |]T3[. The dimensionless lattice Hamiltonian is 

W = -^H = J2 L\n) + ^(-l)"0t(^)0(^) - ^ [0t(n)e'^("V(« + 1) - h-c] (4) 

where a is the lattice spacing, 0(n) is a single-component fermion field at site n, e*^*^"'-' is the 
link variable at link n, and /x = 2m /{g"^ a), x = l/{g'^o?). 

We use a "compact" formulation where the gauge field becomes an angular variable 
< 9{n) < 27r on the lattice, and L{n) is the conjugate spin variable (the lattice equivalent 
of the electric field) 

[9{n),L{m)] = iSnm (5) 

so that L{n) has integer eigenvalues L{n) = 0, ±1, ±2, . . .. 

In the lattice strong-coupling limit x — > 0, the unperturbed ground state |0) has 

L(n) = 0, <j)\nmn) = ^[l-i-in all n (6) 
whose energy we normalize to zero. The lattice version of Gauss' law is then taken as 

L(n) - L{n - 1) = <p\n)<P{n) - ^[1 - (-1)"] (7) 

which means excitations on odd and even sites create =f1 units of flux, corresponding to 
"electron" and "positron" excitations respectively. Eq. (^ determines the electric field 
L{n) entirely, up to an arbitrary additive constant a, which then represents the back- 
ground field. Allowing a to be non-zero, the electrostatic energy term is modified to 
J2n L'^{n) — ^ Y.n{.L{n) + aY . The physics of the background field then matches precisely with 
the continuum discussion. Physics is then periodic in a with period 1, and the background 
field variable is ^ = 27ra . 
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In the weak-coupling limit m/g ^ oo and for background field a = 1/2, there are two 
degenerate vacuum or 'loop' states, ||) and | — |), corresponding to {L{n)+a = ±1/2, all n}, 
respectively. On a finite lattice, the eigenstates will be the symmetric and anti-symmetric 
combinations of those, and we will denote the energy gap between them by Aq. The single- 
fermion state will consist of a single electron with L{n) + a = +1/2 to its left and —1/2 to 
its right; the energy gap between this and the ground state will be denoted Ai. Finally, the 
energy gap to the lowest 2-particle electron-positron state will be denoted A2. 

We will also study two order parameters which can be used to characterize the phase 
transition at 6' = tt. The first one is the average electric field 

r" = ^(E(^H + «))o, (8) 

which in the weak-coupling limit m/g 00 takes values ±1/2 for the zero-loop and one-loop 
states respectively. The second order parameter is the axial fermion density suggested by 
Creutz [| 

= (#7'^/^7)o (9) 



N 



iY^i-ir <P\n)<P{n + l)-h.c. )o. (10) 



Our results are based on the density matrix renormalization group (DMRG) method 
]. The method employed here is the "infinite system" DMRG method, as prescribed by 
White p, used both with open and periodic boundary conditions (OBC and PBC). Due 
to the presence of the electric field on the links, and the differing nature of the odd and 
even-numbered sites, some modifications were made to the method, although nothing that 
changes the spirit of the DMRG. 

For a particular spin configuration with OBC, if we specify the incoming electric field 
for the first site of the chain, then according to Gauss's law the electric field for all the links 
can be deduced. If PBC are imposed, then as there is no particular link to fix the electric 
field, we can have loops of electric flux extending throughout the ring. In the presence of a 
background field we simply add (or subtract) a from the values of the electric fields due to the 
spin configuration. A cutoff, or maximum loop value was chosen such that full convergence 
was reached to machine precision. A loop range of [-5,5] was more than sufficient in most 
cases. 

In the standard "infinite system" method one splits the chain (or ring) into two blocks 
and two sites, and in each DMRG iteration the blocks increase in size by a single site. Due 
to the differing nature of the odd and even sites in the lattice formulation (^, we modify 
this so that two sites are augmented each time, so that the superblock grows by four sites 
in a single DMRG iteration. We truncate the basis in each DMRG iteration such that there 
are a maximum of 125 states per spin sector per loop in a block. 

The accuracy of the DMRG calculation is strongly dependent on the parameter values x 
and m/g. In the worst case scenario {m/g = 0.3; x = 100; 6 = n; N = 256) with PBC, which 
represents the smallest lattice spacing used for this study, and lies in the critical region, the 
ground state energy can be resolved to 1 part in 10^, while the two-particle gap is resolved 
to 1 part in 10^, while the order parameters are accurate to 3 figures. It is more usual to 
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obtain ground state eigenvalues to near machine precision, and "vector" gap errors of 1 part 
in 10^ 

We begin with the case of zero background field, ^ = 0, in order to demonstrate the power 
of the DMRG method. The most accurate results to date are those of Sriganesh et al. ||10|| , 



who used exact diagonalization on lattices up to = 22 sites. With DMRG it is possible 
to go to = 256 sites, such that there is essentially no extrapolation necessary to obtain 
the bulk limit N — > oo. Fig. |^ shows the data for m/g = 0, which is the exactly soluble 
case. For x = 4 through to a; = 100 we have 6 figure convergence. The final continuum limit 
a — > or X — > cxD is obtained by making polynomial fits in powers of 1/a/x = ga as used in 
Ref. [|10|. Our final estimate of / g = 0.56419(4) for this case agrees extremely well with 
the analytic result of = 0.5641896.., and is about 25 times more accurate than the 

corresponding estimate 0.563(1) by Sriganesh et al. [|T0|. 

We now turn to the case of background field 6 = tt. First, we use finite-size scaling 
theory to estimate the position of the critical point, by calculating "pseudo-critical points" 
at each lattice size A^, and each lattice spacing x. The pseudo-critical point is defined [|l^ as 
the point where two successive finite-lattice energy gaps scale like 1/A^. For this exercise we 
used the 'loop gap' Ao/g, which collapses to zero at the pseudo-critical point. The pseudo- 
critical points converge very rapidly to the bulk limit, like 1/A^^. The results are shown 
in Fig. H A quadratic fit in 1/y/x extracts the continuum limit, which we estimate to be 
{m/g)(. = 0.3335(2). This is consistent with the previous estimates by Hamer et al. Q of 
{m/g)c = 0.325(20), or Schiller and Ranft |Tl[, {m/g)c = 0.31(1), but with two orders of 
magnitude improvement in accuracy. 

Finite-size scaling theory |T^ also allows us to estimate the critical indices for the model. 
The order parameters and F^^, for instance, are expected to scale like ~ N~^^'^ [15]. Then 
the logarithmic ratios 

ln[F^((m/^)-v)/F;v_i((m/^)^)] _^ 

\n[N/{N-l)] ^ V ^ ' 

give a direct estimate of the index ratio. Fig. ^ shows the resulting estimates for /9/z/, ob- 
tained from the electric field order parameter, and for 1/z/, obtained from the beta function. 
We see essentially no variation in the exponents with lattice spacing, to within the accuracy 
of our calculations. Our best estimates for the critical exponents are thus z/ = 0.99(1), 
13 /u = 0.125(5). These results provide reasonably conclusive evidence that the Schwinger 
model transition at 6 = ir lies in the same universality class as the one dimensional transverse 
Ising model, or equivalently the 2D Ising model, with {u = 1,13 = 1/8). 

We now turn to estimating continuum limit values for the energy gaps and order param- 
eters. The convergence at ^ = vr is not so good as at 6' = 0. Fig. |^ displays our final results 
for the energy gaps Ao/g between the "loop" states, Ai/g in the 1-particle sector, and A2/g 
in the 2-particle sector, for all values of m/g. 

We see that all gaps vanish at the critical point. The 'loop' gap is zero for all m/g > 
{m/g)c, as predicted by Coleman 0. The 2-particle gap A2/g vanishes at the critical point 
{m/g)c, but on either side of the critical point there is a finite gap, and an almost linear 
behavior with m/g. Note, however, that the behavior is not exactly linear. The 1-particle 
gap vanishes for m/g < {m/g)c, while for m/g > {m/g)c it is very close to half the 2-particle 
gap. Once again, the behavior is very nearly linear in m/g. 
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The pattern of eigenvalues exhibited in Fig. 13 bears an extraordinary resemblance to 



that of the transverse Ising model [|16l. In particular, we see that the energy of the 1- 
particle or 'kink' state vanishes at the critical point, and then remains degenerate with the 
ground state for (m/g) < {m/g)^. Assuming this degeneracy is exact, this indicates that a 
'kink condensate' will form in the ground state for small mass, as discussed by Fradkin and 
Susskind [0. It also indicates the existence of a new 'dual symmetry' in the model, which 
is spontaneously broken in the low-mass region, and has not been explored hitherto. There 
should also be a 'dual order parameter' associated with this symmetry. 

We may also obtain estimates for the order parameters and F° as functions of m/g. 
Our results are displayed in Fig. |^. Both order parameters are zero, within errors, for 
m/g < {m/g)c. It appears that both order parameters turn over and drop abruptly to zero 
as the critical point is approached from above, consistent with the small exponent u = 1/8 
found previously. The axial density F^ decreases steadily towards zero at large m/g, whereas 
F" approaches the expected value of 1/2. More detailed results will be given in a forthcoming 
paper . 
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FIGURES 
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FIG. 1. A configuration of "half-asymptotic" charged fermions at background field F = g/2. 
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FIG. 2. "Vector" mass gaps for m/g = 0, finite lattices N = 10-256, and various lattice spacings. 
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FIG. 3. Critical line in the m/g versus 1 / \/x = ga plane. Open circles are our present estimates, 
and squares are the previous results of Hamer et al. which are in good agreement. The dashed 
line is a quadratic fit to the data in ga. 
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FIG. 4. 'Logarithmic' ratio estimates of critical indices —1/u and — /J/i' for lattice spacing 
l/\/x = ga = 0.45. Quadratic fits in 1/N provide the bulk extrapolations. We estimate here 
= 1.00(2) and (3/i^ = 0.125(5). 
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FIG. 5. Final estimates for gaps in the 0-particle, 1-particle and 2-particle sectors at = tt. 
Dashed lines are merely to guide the eye. 
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FIG. 6. Order parameters = ((L + a))o , F^ = {iil^^5il^/g)o near the critical region. Dashed 
hnes are merely to guide the eye. 
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